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ABSTRACT 

We report results for the angular three-point galaxy correlation function in the APM survey 
and compare them with theoretical expectations. For the first time, these measurements extend 
to sufficiently large scales to probe the weakly non-linear regime. On large scales, the results 
are in good agreement with the predictions of non-linear cosmological perturbation theory, for a 
model with initially Gaussian fluctuations and linear power spectrum P(k) consistent with that 
inferred from the APM survey. These results reinforce the conclusion that large-scale structure 
is driven by non-linear gravitational instability and that APM galaxies are relatively unbiased 
tracers of the mass on large scales; they also provide stringent constraints upon models with 
non-Gaussian initial conditions and strongly exclude the standard cold dark matter model. 



Subject headings: cosmology: theory; cosmology: observations; cosmology:large-scale structure 
of universe; methods: numerical; methods: analytical 
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1. 



Introduction 



Traditionally, two-point statistics, the auto-correlation function £(r) and the power spectrum P(k), 
have been the dominant benchmarks for testing theories of structure formation. However, with the advent of 
large galaxy surveys and the development of non-linear cosmological perturbation theory and large N-body 
simulations, it has become clear that higher-order correlations provide new probes of large-scale structure. 
In particular, the N > 2-point functions and moments test models for bias — the relation between the galaxy 
and mass distributions — and constrain non-Gaussianity in the initial conditions (Fry & Gaztahaga 1993, 
Frieman & Gaztahaga 1994, Gaztahaga 1994, Gaztahaga & Frieman 1994, Fry 1994, Fry & Scherrer 1994, 
Jaffe 1994, Fry, Melott, & Shandarin 1995, Juszkiewicz, et al. 1995, Chodorowski & Bouchet 1996, Gaztahaga 
& Mahonen 1996, Szapudi, Meiksin, & Nichol 1996, Jing 1997, Matarrese et al. 1997, Mo, Jing, & White 
1997, Verde et al. 1998, Scoccimarro et al. 1998, Scoccimarro, Couchman, & Frieman 1998). 

The galaxy three-point function £ has been measured in several angular and redshift catalogs (Peebles 
& Groth 1975, Peebles 1975, Groth & Peebles 1977, Fry and Scldner 1982, Bean et al. 1983, Efstathiou 
& Jedredjewski 1984, Hale-Sutton ct al. 1989, Baumgart & Fry 1991, Jing, Mo, & Bocrncr 1991, Jing & 
Boerner 1998). Given the limited volumes covered by these surveys, these measurements generally probed the 
three-point function on scales 0.1 < r < 10 h^ 1 Mpc. They established that galaxies cluster hierarchically: 
on these scales, the hierarchical three-point amplitude, defined by 



(with Xij = \xi — Xj\) is nearly constant, independent of scale and configuration, with values in the range 
Q ~ 0.6 — 1.3, depending on the catalog. This hierarchical form is consistent with expectations from N- 
body simulations in the non-linear regime (e.g., Matsubara & Suto 1994, Fry, Melott, & Shandarin 1993, 
Scoccimarro et al. 1998, Scoccimarro & Frieman 1998). 

On larger scales, r > 10 Mpc, in the weakly non- linear regime where £(r) < 1, non- linear pertur- 
bation theory (PT) — corroborated by N-body simulations — predicts that Q becomes strongly dependent on 
the shape of the triangle defined by the three points x t (Fry 1984, Fry, Melott, & Shandarin 1993, Jing & 
Boerner 1997, Scoccimarro 1997, Gaztahaga & Bernardeau 1998, Scoccimarro et al. 1998). The large-scale 
configuration-dependence of the three-point function is characteristic of the non-linear dynamics of gravita- 
tional instability and is sensitive to the initial power spectrum, to the bias, and to non-Gaussianity in the 
initial conditions. Thus measurement of the three-point function and its shape dependence should provide 
a powerful probe of structure formation models. Measurements of higher-order galaxy correlations on large 
scales have been confined so far to volume-averaged correlation functions — the one-point cumulants Sn — 
and their cousins (e.g., cumulant correlators). These statistics, based on counts-in-cells, are computationally 
easier to measure than the corresponding N— point functions and, being averages over many N— point con- 
figurations, can be measured with higher signal to noise as well. On the other hand, this averaging process 
by definition destroys the information about configuration-dependence contained in the TV— point functions 
themselves. 

In this Letter, we report measurement of the projected three-point function on large scales in the APM 
Galaxy Survey (Maddox et al. 1990), compare the results with predictions of non-linear perturbation theory 
(PT) and N-body simulations, and briefly discuss their implications for bias and non-Gaussianity. 

We first recall the perturbative expression for the three-point function. To lowest order, the statistical 
properties of the density contrast field, S(x) = (Sp(x)/p) — 1, are characterized by the auto-correlation 
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function £,(xi 2 ) = (8(xi)8(x 2 )), and its Fourier transform, the power spectrum P(k), where (8(ki)8(k 2 )) = 
(2Tr) 3 P(k 1 )Sl(k 1 + fc 2 ); explicitly, 

^(r) = ^ J dkk^P(k)^l . (2) 

In the results to be shown below, we consider the cold dark matter (CDM) family of models, with shape 
parameter T = £lh ranging from 0.5 (standard CDM) to 0.2. We also consider a model consistent with the 
linear power spectrum inferred from the APM survey itself (Baugh & Gaztanaga 1996), 

PAPM - likc(fc) = [i + (ktomrr ' (3) 

with A: in ft, Mpc -1 . For Gaussian initial conditions, to linear order in 8, the connected N > 2-point 
functions vanish. In second-order perturbation theory, the three-point correlation function C( a3 ij x 2, ^3) = 
(8(xi)8(x 2 )8(x 3 )) can be expressed as (Fry 1984, Jing & Boerner 1997, Gaztanaga & Bernardeau 1998) 
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where 



*(r)=^ ldkP(k)°-^l, (5) 
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and/'(x) = df/dx. In eqns.(^)-(j^), the power spectrum and auto-correlation function are implicitly evaluated 
at linear order in perturbation theory; we expect the leading-order result (Q) to be valid in the weakly non- 
linear regime (£ <C 1) but to break down for £(r) > 0.5. Eqn.Q) was derived for the Einstein-de Sitter 
(fi m = 1) model, but it is known to be an excellent approximation for fl > 0.1 (Bouchet et al. 1992, 1995, 
Bernardeau 1994, Catelan et al. 1995, Scoccimarro et al. 1998, Kamionkowski & Buchalter 1998). For non- 
Gaussian initial conditions, there will in general be corrections to (Q) which depend on the initial (linear) 
three- and four-point functions (Fry & Scherrer 1994). 

The expressions above apply to unbiased tracers of the density field; since galaxies of different types 
(e.g., ellipticals and spirals) have different clustering properties, we know that at least some galaxy species 
are biased. As an example, suppose the probability of forming a luminous galaxy depends only on the 
underlying mean density field in its immediate vicinity. Under this simplifying assumption, the relation 
between the galaxy density field 8 g (x) and the mass density field 8(x) is 8 g (x) — f(8(x)) = Yl n b n 8 n ) where 
b n are the bias parameters. To leading order in perturbation theory, this local bias scheme implies 

Ux) = bi^x) , (6) 

C fl (sBi,X2,a!3) = b\Q(x 1 ,x 2 ,x 3 ) + b\b 2 [t(x l2 )S,(x 13 ) + £(x 12 )£(x 23 ) + £(x 13 )£(x 23 )} , (7) 
and therefore (Fry & Gaztanaga 1993, Fry 1994) 
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Gaztanaga & Frieman (1994) have used the corresponding relation for the skewness S3 to infer 61 ~ 1, 62 — 
from the APM catalog, but the results are degenerate due to the relative scale- independence of 53. [j] On 
the other hand, Fry (1994) used the projected bispectrum from the Lick catalog to infer b\ ~ 3; however, in 
order to extract a statistically significant Q g from the catalog, an average over scales including those beyond 
the weakly non-linear regime was required. As we will see below, the configuration-dependence of Q g on 
large scales in the APM catalog is quite close to that expected in PT, suggesting that b\ is of order unity 
for these galaxies. The simple model above undoubtedly does not capture the full complexity of biasing 
(e.g., Mo & White 1996, Blanton et al. 1998, Dekel & Lahav 1998, Sheth & Lemson 1998), but it provides 
a convenient framework that is well matched to the quality of the current data. 

In a projected catalog with radial selection function <j){x) (normalized such that J dxx 2 <f>(x) = 1), the 
galaxy angular two- and three-point functions at small angular separations (9ij <C 1) are given by (Peebles 
& Groth 1975, Peebles 1980) 

p OO pOO 

w g (d) = 2 dxx i <p 2 (x)F- 2 (x) du£ g {r;z) (9) 
Jo Jo 

/•oo />oo />oc 

z g (Oi2,O 13 ,023)= dxx 6 cf> 3 (x)F- 3 (x) / dudv C g (n 2 , ria, r 23 ; z) , (10) 

J0 J-00 J-00 

where #12, #13, and #23 are the sides of a triangle projected on the sky. We assume the proper separations 
rij are small compared to the mean depth of the sample; in this case, 

n^^lik/Ff + xHl] 112 , (ii) 

with /x2 = u, fi3 = v, and /2s = u — v. Here, x is the comoving radial (coordinate) distance to redshift 
z, and F(x,SY) is a geometrical factor which relates proper and coordinate distance intervals, F(x) = [1 — 
(Hox/c) 2 ^-!)} 1 / 2 . 

The radial selection function for the APM Galaxy Survey can be approximated by (Gaztahaga & Baugh 
1998) 

cf>(x) AP M =Cx- b ex V ~ x2 /° 2 , (12) 

with b = 0.1 and D = 335/i _1 Mpc. At this depth, to an accuracy of better than a few percent we can 
approximate F by its Einstein-de Sitter value F = 1. The projected hierarchical amplitude is defined by 
analogy with eqn. (||) , 

„ (a a a \ — z g (9 12 , gl3j ^23) 

93^12, 013, V23) = — ttp^; j- a — r- , a \ 77^- ■ (4 a J 

W g [B\2 )Wg [Vl3 ) + W g [Oil )Wg (023 ) + W g (#i 3 )w g (#23 ) 

In projecting eqns.(j^) and (^) in PT we assume leading-order perturbative growth for the redshift-evolution 
of £ and £ for f2 m = 1; in this case, both Q and 53 are independent of the power spectrum normalization 
(e.g., independent of erg). 

For the analysis, we use the equal-area projection pixel map of the APM survey, with an area roughly 
120° x 60°, containing N ga i ~ 1.3 x 10 6 galaxies brighter than bj = 20 and fainter than bj = 17. Each pixel 
has an area (O p ) 2 — (0.06) 2 sq. deg., and the mean galaxy count per pixel is (N) = N gal /N pix = 0.97. The 



1 Bernardeau (1995) pointed out a systematic correction to the way the Sj predictions should be projected. After taking 
into account the correct selection function and the uncertainties in the APM shape of P(k), this effect is less significant that 
claimed by Bernardeau, and the original interpretation is still valid, in a greement with the results for 53 presented here. 
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estimator of the density fluctuation amplitude at the ith pixel is Si = (Ni/(N)) — 1, where 2Vj is the galaxy 
count in that pixel. The estimator for the galaxy angular two-point function is then (Peebles & Groth 1975) 

^) = ^I>Wi(0). (14) 

i,j 

where Ng = . Wij (0) is the number of pairs of pixels at separation 9 in the survey region, and the 
angular window function Wij (9) = 1 if pixels i and j are separated by \9i — 8j\ = 9 ± © p , and otherwise. 
The reduced angular three-point function is estimated as 

2(012,013,023) = -^J2Si6 j 5 k W m (0 12) e 1 s,e23) , (15) 

where Ng — - k Wijk is the count of triplets of pixels with angular separations 0i2, #13, 023 in the survey 
region, and Wijk = 1 f° r pixel triplets with that angular configuration and otherwise. We note that, in 
the limit of counts of pairs and triplets of objects, these estimators are equivalent to the minimum variance 
estimators of Landy & Szalay (1993) and Szapudi & Szalay (1998). We employ these estimators for angular 
separations > 0.5 deg, at least an order of magnitude larger than the pixel scale; in this limit, finite pixel-size 
corrections should be less than a few percent and have been neglected. 

To test the validity of PT on the scales of interest, to verify the algorithm for measuring the angular 
three-point function, and to check projection and finite sampling (and boundary) effects in the APM survey, 
we use the simulated APM maps of Gaztahaga & Baugh 1998. These are created from N-body simulations of 
the SCDM (with <j 8 = 1) and APM-like (cr 8 = 0.8) models, with box size of 600 h^ 1 Mpc and 128 3 particles, 
both with O m = 1. From each N-body realization we make 5 APM-like maps from 5 different observers; 
the observers are spaced sufficiently far apart compared to D that the 'galaxies' they observe do not appear 
to be strongly correlated (Gaztahaga & Bernardeau 1998). As the simulation is done in a periodic box, we 
replicate the box to cover the full radial extent of the APM (around 1800 hr 1 Mpc, at which distance the 
expected number of galaxies is smaller than unity). To account for possible boundary effects, we employ 
the APM angular survey mask, including plate shapes and holes. We have also made larger angular maps 
without the APM mask and find no significant differences from the results shown below, from which we 
conclude that the APM mask does not affect the estimations of the 2 or 3-point functions on the scales 
under study here. The projected simulations have ~ 30% lower mean surface density of galaxies than the 
APM. We have diluted the APM pixel maps by this amount (by randomly sampling the galaxies) and found 
that this dilution has a negligible effect on the clustering properties within the errors. 

To estimate the projected 3-point function we need to find all triplets in the pixel maps which have a 
given triangle configuration and repeat the process for all configurations. Consider a triplet of pixels with 
labels 1,2,3 on the sky. Let 0i2 and #13 be the angular separations between the corresponding pairs of pixels 
and a the interior angle between these two triangle sides. One can characterize the configuration-dependence 
of the three-point function by studying the behavior of 53 (a) for fixed 0i2 and $13 . Our algorithm for counting 
such triplets is as follows. For each pixel in the map, we find all pixels that lie in two concentric annuli of 
radii #12 ± 0.5 and #13 ± 0.5 about it, where 9ij is now measured in units of O p . We count all pairs of pixels 
between the two annuli; this requires (m\i x mx3)/2 operations, where rriij — 2ir9ij. As seen from the central 
pixel, each pair has an angular separation a, and the results are binned in a. We repeat this procedure for 
each pixel in the map, building the estimators ( |l4| ) and ([[5]). Thus, to find all triplets with separations #12 
and 0i3 requires only N pix (mi2 x mxz)/2 operations; for 9 12 = 0i3 = 1 deg, this computation only takes a 
few hours of CPU time on a modest workstation. By contrast, a naive 0{N^ ix ) operation with N pix ~ 10 6 
would be quite a lengthy computational task given current computer power. 
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The two lower panels in Fig. |T| show results for 93(a) for the SCDM and APM-like models, for 812 = 
6*i3 = 2 degrees, projected at the depth of the APM survey. The configuration-dependence of the hierarchical 
three-point amplitude is seen to be quite sensitive to the shape of the power spectrum. Both the shape and 
amplitude of (73(a) predicted by PT (solid curves) are reproduced by the N-body results (points) even on 
these moderately small scales (at the mean depth of the APM, 2 deg corresponds to ~ 14 h _1 Mpc). Part of 
this agreement traces to the fact that (73 involves a weighted sum over spatial 3-point configurations covering 
a range of scales; due to the shape of the power spectrum, Q increases on large scales, so these configurations 
(which are further in the perturbative regime) are more heavily weighted in projection. The error bars 
on the simulation results are estimated from the variance between 10 maps (5 observers each in 2 N-body 
realizations), assuming they are independent, and correspond to the 1-a interval of confidence for a single 
observer (i.e., they are not divided by — 1). 

The top panel of Fig. [I] and all those of Fig. || show the measurements of (73(a) in the APM survey itself 
at 812 = #i3 = 0.5 — 4.5 degrees. Closed squares correspond to estimations in the full APM map. The values 
at a = are in agreement with the cumulant correlators, C12 = 9(0), estimated (with 4x4 bigger pixels) by 
Szapudi & Szalay (1999). The mean values are comparable to the values of S3/3 in the APM (Gaztanaga 
1994, Szapudi et al. 1995). Also notice that at the scales considered here, > 0.5 deg, the values of S3 in the 
APM and the EDSGC (Szapudi, Meiksin, & Nichol 1996) are very similar (Szapudi & Gaztanaga 1998). 

The APM results are compared with the values of (73 for the APM-like spectrum in PT (solid curves) 
and in simulations (open triangles with errorbars). As the APM-like model has, by construction, the same 
w(6) as the real APM map, we assume that the sampling errors should be similar in the APM and in the 
simulations. This might not be true on the largest scales, however, where systematics in both the APM survey 
and the simulations are more important (e.g., the simulation might be affected by the periodic boundaries, 
and the power in the simulation and in the survey may differ on the largest scales). At scales 8 > 1 deg, 
the agreement between the APM-like model and the APM survey is quite good; this corresponds roughly 
to physical scales r > 7 h^ 1 Mpc, not far from the non-linear scale (where £2 — 1)- At 4.5 deg, the APM 
agrees better with the PT prediction than with the simulation results, which show large variance; this could 
be an indication that the 600 h _1 Mpc simulation box is not big enough, leading to larger sampling errors 
than the APM itself. We also note that the SCDM model clearly disagrees with the APM data for 93; this 
can be seen at 9 12 = 6*13 = 2 deg (Fig.|l|) and for 6*12 = 6*13 = 3.5 deg (dashed curve in lower-left panel of Fig. 
^). This conclusion is independent of the power spectrum normalization. 

At smaller angles, 8 < 1 deg, 93 in the simulations is larger than in either the real APM or PT (top-left 
panel in Fig. |^). The discrepancy between simulations and PT on these relatively small scales is clearly 
due to non-linear evolution. The interpretation of the discrepancy with the real APM is less clear: a 
number of assumptions underlying the simulations could affect the final results at these non-linear scales. 
For example, systematic uncertainties in the APM selection function or a linear bias would lead one to infer 
a different linear APM-like power spectrum from the w(6) data. Also, a model with low f2 m would undergo 
less non- linear evolution, which might give a better match to the APM results for 93; this could provide an 
interesting test for a low-density universe. This discrepancy at non-linear scales is similar to the one found 
in Baugh & Gaztanaga (1996), where the real APM values of Sj were closer to the PT predictions than 
to simulations. Other possible contributions to this effect include non-linear bias and non-linear projection 
effects (Gaztanaga & Bernardeau 1998). 

The open circles in each figure show the mean of the estimations of 93 in 4 disjoint subsamples of the 
APM survey (equally spaced in right ascension, as in Baugh & Efstathiou 1994). For illustration, the values 
of (73 for each of the 4 zones are shown in the top panel of Figure [|: the dotted, short-dashed, long-dashed, 
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and dot-dashed curves correspond to zones of increasing RA (the middle two of these correspond to relatively 
lower galactic lattitude). These estimations of (73 are subject to larger finite- volume effects, because each 
zone is only 1/4 the size of the full APM. Because the zones cover a range of galactic latitude, a number of 
the systematic errors in the APM catalog (star-galaxy separation, obscuration by the galaxy, plate matching 
errors) might be expected to vary from zone to zone. We find no evidence for such systematic variation 
in (73: the individual zone values are compatible with the full survey within the (sampling) errors in the 
simulations (compare the top and middle panels in Figure [l]) . On larger scales, 9 > 3 deg, the individual 
zone amplitudes exhibit large variance, and boundary effects come into play. 

On all scales, we find large covariance between the errors of 93: the data points at different a are 
strongly correlated. This is illustrated in the middle panel of Fig. |l]: the dotted and continuous curves 
correspond to results for 2 of the 10 observers. Sampling or finite-volume effects are seen to produce a 
systematic vertical shift in the curve rather than a scatter around some mean value. A similar trend is 
found in the real APM catalog (compare the 4 curves in the top panel of Fig. |l|). This covariance can be 
studied analytically (Hui, Scoccimarro, Frieman, & Gaztanaga, in preparation) and comes from fluctuations 
on scales comparable to the sample size. Ratio bias and integral-constraint bias (e.g., Hui & Gaztanaga 
1999) could also be important. 

Other possible sources of systematic discrepancies between the model predictions and the APM results 
include the shape of the APM radial selection function, the evolution of clustering, and the shape of the linear 
power spectrum. We find that the first two effects introduce differences smaller than 10% in the amplitude 
of (?3 (in agreement with Gaztanaga 1995), which are not significant given the errors. The uncertainty in 
the shape of the linear P{k) is more important, and as mentioned above is critical for the interpretation of 
(73 at the smallest angles. Nevertheless, at large scales (9 > 1 deg) these uncertainties appear to be within 
the errors (when we take into account that the errors are strongly correlated). This is illustrated in the 
bottom-left panel of Fig. || — the two solid curves bracketing the APM-like results show the PT predictions 
for two power spectra which conservatively bracket the uncertainties in the linear spectrum inferred from 
the APM w{9) (see Fig.13 of Gaztanaga & Baugh 1998): P(k) ~ k a /[l + (fc/0.06) 3 ], with a = 0.2 and 1.2. 
The model with a = 1.2 (the lower solid curve at small a) appears to give a better match to the g3 results 
in the APM than the central APM-like model of Eq.||; for reference, a CDM model with T = 0.3 gives a 
nearly identical value for (73(a) on these scales. Thus, although the APM results for q% generally fall below 
the PT predictions on angular scales 9 > 2 deg, they are consistent within the sampling errors and given the 
uncertainties in the shape of the linear power spectrum inferred from the APM w(9). 

To place accurate constraints upon bias models and initially non-Gaussian fluctuations, we must quanti- 
tatively model the covariance between the estimates of (73; this will be done elsewhere, but we can nevertheless 
get a qualitative sense of the limits here. We expect the strongest constraints to come from intermediate 
scales, 9 ~ 1 — 2.5 deg, where both the sampling errors and the non-linearities are small. The upper right 
panel of Fig. || shows the PT predictions for the APM-like model with linear bias parameter b\ = 2 (dashed 
curve) and a non-linear bias model with 61 = 1, 62 = —0.5 (dotted curve). Even if the errors are 100% 
correlated, these models are clearly ruled out by the APM data; we conservatively conclude that b\ < 1.5 is 
required for a simple linear bias model to fit the APM data. Note that a model with bi = 1.5 and b 2 = 0.5 
would have roughly the correct amplitude for q 3 , but its shape would be flatter than the data, especially at 
larger 9. 

As a simple example of a non-Gaussian model, the dotted curve in the lower-left panel of Fig. || 
shows the leading-order prediction for the x 2 isocurvature model (Peebles 1997, 1998a, 1998b, Antoniadis, 
ctal. 1997, Linde & Mukhanov 1997, White 1998) with the APM-like spectrum. In this model, the initial 
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density field is the square of a Gaussian random field, and the leading-order 3-point function is simply 
C = 2[2^(a;i2)^(xi3)^(a;23)] 1 ^ 2 . Clearly, the projected three-point function for this model is substantially 
larger than that of the corresponding Gaussian model for intermediate a; both the amplitude and shape are 
discrepant with the APM data. To make this comparison precise, the non-linear corrections for this model 
should be self-consistently included; however, these corrections are expected to increase rather than reduce 
the 93 amplitude, likely making the disagreement worse. In general, if we assume no biasing, the initial 
non-gaussianities are restricted to | (73 1 < 0.5. 

At angular scales 9 > 1 deg, corresponding to physical scales for which £ < 1, the agreement between 
PT and the APM survey for the angular three-point amplitude (73 is quite good, implying that APM galaxies 
are not significantly biased on these scales and that their spatial distribution is consistent with non-linear 
evolution from Gaussian initial conditions. This substantiates and extends the conclusions of Gaztanaga 
(1994,1995) and Gaztanaga & Frieman (1994). 
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Fig. 1. — The two lower panels show perturbation theory predictions (solid curves) and N-body results (points 
with sampling errors) for the projected 3-point amplitude 93(a) at fixed 812 = #13 = 2 deg for a survey with 
the APM selection function. N-body results are the mean from 10 observers. Middle panel corresponds to 
the initial APM-like power spectrum of Eqn.(H) and lower panel to the SCDM model. Additional curves in 
the middle panel show results for two different observers. The top panel compares the PT prediction for the 
APM-like model with the real APM measurements (closed squares and open circles correspond to the full 
APM map and to the mean of 4 zones; other curves in the top panel show results for each of the zones). 




Fig. 2. — The projected 3-point amplitude qz in PT (solid curves) and N-body results (open triangles with 
errorbars) for the APM-like power spectrum are compared with (73 measured in the APM survey (closed 
squares and open circles, with same meanings as in Fig. 1). Each panel shows the amplitude at different 
#12 = 013- In upper right panel, dotted and dashed curves correspond to PT predictions with 61 = 1,6 2 = 
—0.5 and b\ = 2, 62 = 0, respectively. In the lower left panel, upper and lower solid curves conservatively 
bracket the uncertainties in the inferred APM-like power spectrum, long-dashed curve corresponds to SCDM, 
and the dotted curve shows the leading-order prediction for the x 2 non-Gaussian model. 



